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Abstract 

The Sivers parton distribution function has been predicted to obey a particular "universahty relation" , 
namely to have opposite sign in semi-inclusive deeply inelastic scattering (SIDIS) and the Drell-Yan 
process. We discuss how, on the basis of present HERMES data, this remarkable prediction of the QCD 
factorization approach to the description of single spin asymmetries related to the Sivers effect could be 
checked experimentally in future experiments at PAX and COMPASS. 

1 Introduction 

It was understood early [1] that single spin asymmetries (SSA) in hard processes, such as those observed 
in p^p ttX [2, 3] or in SIDIS [4, 5, 6, 7, 8, 9], cannot be explained by means of leading twist coUinear 
QCD factorization. One of the non-perturbative effects which could account for such SSA considers a non- 
trivial correlation between (the transverse component of) the nucleon spin St and intrinsic transverse parton 
momenta px in the nucleon [10], and is quantified in terms of the so-called Sivers function f^{x,p'^) [11]. 
The effect is referred to as "(naively) T-odd", since it is proportional, e.g., in the infinite momentum frame 
where the nucleon momentum P^r — > oo, to the T-odd structure (St x pT)PAr. The Sivers effect was shown 
to be able to explain the SSA in p^p ttX [12], though also other mechanisms exist which could contribute 
in this reaction [13, 14, 15]. 

The precise definition of f^{x,Px) in QCD was worked out only recently [16, 17, 18]. A particularly 
interesting feature of the Sivers function concerns its universality property. This property ensures for usual 
parton distributions that one deals with, e.g., the same unpolarized parton distribution fi{x) in SIDIS 
and in the DrcU-Yan process (DY): /i(a;)siDis — fi{x)DY- In tiie case of the Sivers function (and other 
T-odd distributions) the universality property takes, however, a different form. On the basis of time-reversal 
arguments it was predicted [17] that in SIDIS and DY have opposite sign, 

/it(^i Pt)sidis = ^/it(^j Pt)dy • (1) 

The experimental check of Eq. (1) would provide a thorough test of our understanding of the Sivers effect 
within QCD and, hence, our understanding of SSA. It would crucially test the factorization approach to the 
description of processes sensitive to transverse parton momenta [19, 20, 21]. 

In this work we shall discuss how the relation (1) could be checked experimentally in the Polarized 
Antiproton eXperiment (PAX) planned at GSI [22, 23]. A primary goal of this experiment will be to 
provide a polarized antiproton beam and to measure the transversity distribution h\{x)^ c.f. [24]. However, 
PAX will also be well suited to access the Sivers function via SSA in pp^ — > ji^ ji" X or p^p pt^^'X 
[22, 23]. In the COMPASS experiment at CERN [25], making use of a tt~ beam, one would also be able to 
study the Sivers function via SSA in tt~p^ pt^ pL~X. 

In order to estimate the magnitude of the Sivers effect in those experiments we will roughly parameterize 
f^{x, Py)siDis from the (preliminary) HERMES data [7] using as a guideline relations derived from the QCD 
limit of a large number of colours Nc [26]. Such large- iVc relations are observed to hold in nature within 
their expected accuracy [27] and, as a byproduct of our study, we shall observe that this is also the case 
here. On the basis of the obtained parameterization we estimate SSA in the Drell-Yan process for the PAX 
and COMPASS experiments. We also comment briefly on parameterizations of reported previously in 
the literature and on model calculations. 



2 The Sivers function 



A definition of the unintegrated unpolarized distribution function fi{x, Py) and the Sivers function f^{x, p^) 
can be given in terms of the light-cone correlator 



$3(x,pt) = y^^'^e'f«(P,5T|V^q(0)7^n'i W[0,$; process] ^,(0|P,5t) 
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= mx,p'T)+ftT'i^,P'T) , (2) 

where the dimensionless light-like vectors n± are defined such that n+ • n_ = 1. (See Ref. [28] for a precise 
definition and the meaning of unintegrated distribution functions in QCD.) 

The Wilson link W[0, ^; process] is defined in Fig. 1, c.f. Refs. [17, 18]. For observables integrated over 
Pt the process dependence of the gauge link usually cancels out. However, the situation is different for f^. 
If one neglected the gauge link, under time-reversal the Sivers function would transform into its negative, 
i.e. it would vanish [14]. However, initial or final state interactions [16, 29], needed to obtain non-zero SSA 
[30], generate a Wilson link for the Sivers function in any gauge [17, 18]. Under time reversal the gauge link 
of SIDIS is transformed into the gauge link of DY. and vice versa [17]. (The gauge link structure in the more 
involved hadronic process p^p nX was addressed in [31].) This yields the peculiar universality relation in 
Eq. (1). 

Only little is known about the non-perturbative properties of the Sivers function. In Ref. [32] bounds 
were derived from the positivity of the spin density matrix, which constrain in terms of other transverse 
momentum dependent distributions including also so far experimentally unknown functions. Eliminating 
the unknown distribution, at the price of relaxing the bound, one obtains 

^lAV(^,pi.)|</f(a;,pi.). (3) 

The average parton transverse momentum defined as (p^) = Jdx J d^pT Pt^'^{x, Pt) was shown [33, 34] 
(with analogously defined gluon transverse momentum) to obey the relation 

E (Pr)=0. (4) 

a=g,u,d, ... 

Eq. (4) is more than trivial momentum conservation in the plane transverse to the hard momentum flow, 

since it connects transverse momenta due to final state interactions which the scattered quark experiences in 
incoherently summed scattering events [33, 34]. Inserting the definition of $^(a;,pT) into Eq. (4) and using 
the fact that (analogously for gluons) 

Jdx J d^pT ^ h'r'ix, p't) = S'' J dx ft-^^^ix) , with A-^^(^)«(x) = J d^pT ^ A-^t'(x, p|) , (5) 
one obtains 

E [dxf^^'^^{x)=0. (6) 

a—g,u,d, ... 

The sum rule (6) may prove a useful constraint for parameterizations of the Sivers function, in particular, 
because it is f^^\x) which enters in a model independent way cross sections properly weighted with 

transverse momenta [11]. 

Another property of the Sivers function, which will be used later on, is the relation derived in the limit 
of a large number of colours Nc in QCD [26], namely 

/]^(x,p|n) = — /j^(a;, p|n) modulo corrections. (7) 

It should be noted that in the large- iVc limit xN,. = 0{1), such that Eq. (7) can be expected to be satisfied 
in the valence region of not too small and not too large x to within an accuracy of 0{1/Nc) [27]. 

Neglecting strange and heavier quarks, which is a reasonable assumption in the case of the nucleon, one 
obtains from Eqs. (6, 7) that the Sivers gluon distribution is suppressed in the large-A^c limit with respect 
to the quark distribution functions. More precisely, it is of the same order of magnitude as the fiavour 
singlet combination^, i.e., (/f^ + fir ) ^ fir ^ C'(iV^). Thus, in the large- A^c limit the gluon Sivers effect 

^What matters in the large A'c-counting is the spin-flavour symmetry of the involved operator. In this respect the operator 
entering the flavour singlet and the gluon Sivers function have the same behaviour and thus the same large- AT^. counting. We 
thank Pavel Pobylitsa for discussions on this point. 
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Figure 1: The path of the process-dependent gauge hnk W[0,^; process] which enters the definition of the Sivers function in 
SIDIS and DY. 



can be expected to be suppressed with respect to the non-singlet quark Sivers effect at not too small x, 
which is an interesting constraint for phenomenological studies. In order to obtain a feeling to which extent 
such large-iVc relations may be expected to hold, it is interesting to mention that the hclicity distribution 
fimction exhibits a similar behaviour in the large-A^c limit, namely \{gi — gi){x)\ ~ 0{N^) is larger than 
I (ffi + ffi ) {^) I ^ Iff? i^) I ~ C>{Nc). This is - for quarks - roughly consistent with phenomenology, and predicts 
a suppression of the (presently poorly known) helicity gluon distribution with respect to the unpolarized 
gluon distribution function gf{x)/ff{x) ~ at moderate values of x [27]. 

The Sivers function was studied in several models, in which the gauge link was modeled explicitly by 
considering the perturbative effect of one-gluon exchange. Calculations based on the spectator model [16, 
35, 36] and the bag model [37] yield a sizeable but a negligible and, thus, for the chosen parameter 
sets, do not respect the large Nc counting rule in Eq. (7). 

In a large class of chiral models, which are based on Goldstone boson and (effective) quark degrees of 
freedom, T-odd distributions are zero [38]. This can be understood by recalling that in such models there 
are no gluons, whose presence is crucial to generate T-odd effects via the gauge link structure. Combining 
the no-go-theorem (concerning modeling of in chiral models) of Ref. [38] with notions from the instanton 
model of the QCD vacuum [39] one is lead to the suspicion [40] that T-odd distributions could be suppressed 
with respect to T-even distributions in the instanton vacuum model, which is supported by estimates [41]. 
For a discussion of possible instanton effects in the Drell-Yan process see [42] and references therein. 

By assuming that the SSA in p^p wX [2] is dominated by the Sivers effect phenomenological pa- 
rameterizations of the Sivers function were obtained [12, 43, 44] which, worthwhile stressing, approximately 
respect the large- A^^ pattern in Eq. (7). The assumption that this proc;ess is dominated by the Sivers effect 
seems to be reasonable in the light of recent studies [45, 46], which show that the contribution of the Collins 
effect in this process is small. However, one also has to keep in mind twist-3 effects [13, 15] which might be 
equally important. 

Let us finally mention that a relation of /{^ to the generalized parton distribution E'^{x,^,t) was pro- 
posed, namely the leading light-cone Fock component of the Sivers function may be represented as a convo- 
lution of the Wilson link and the same overlap integrals between light-cone wave functions differing by one 
unit of orbital angular momentum, which enter the description of E'^{x,^,t) [16, 47]. This has been seen 
explicitly in a quark-diquark model calculation [48] and is compatible with the large Ac-hmit in the sense 
that and /{^ have the same large- A^c behaviour [49]. From these relations it was concluded that 

/d^/lT^s!ors(^)<0, Jdxf^^2Li^)>0. (8) 

The above connection to the generalized parton distribution could further be exploited to draw conclusions 
on the large- a; behaviour of the Sivers function. Since E'^{x, ^, t) oc (1 — x)^ at large x [50] one may conclude 
that also f[^'^''{x) oc (1 — x)^. One must keep in mind, however, that the above assumes dimensional 
counting behaviour of the usual quark distribution functions, which yields /'(a;) oc (1 — a;)^ at large x [51]. 
While for /"(x) parameterizations, e.g. [52, 53], are roughly compatible with this prediction, this is not the 
case for ff{x). However, if this were true also for E''{x,^,t) and f[^'^''{x), there need not to be a conflict 
with large- Ac relations, such as Eq. (7), which apply only as long as xN^. = 0{1). 
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Figure 2: Kinematics of the SIDIS process Ip — * I'hX (left), and the Drell-Yan process p^h ^ I'^l X (right) in the lab frame. 

3 Sivers effect in SIDIS 

Consider the process Ip'^ I'hX (see Fig. 2) where denotes the transverse (with respect to the beam) 
target polarization. Let P, I (/') and Ph denote respectively the momentum of the target proton, incoming 
(outgoing) lepton and the produced hadron. The relevant kinematic variables are s := {P + l)^ , q = I — I' 
with = — <7^, X — 1^^, y — z ~ ^^p^- The azimuthal SSA asymmetry is defined as 

iEa^K0.>S.) + A^K'^.,0s. + vr)} ' 

where N^^^^{(j)i, 4>s.i) are the event counts for the respective target polarization (corrected for depolarization 
effects). The z-axis is chosen in direction of the virtual photon (see Fig. 2), in agreement with the HERMES 
convention [4] . The angles iph and ips are the azimuthal angles of the produced hadron and the target spin. 
Neglecting power suppressed terms cx Mj^/Q"^ the SSA (9) is given by [11]^ 



A-^t'-'^ix) = 2 ^-^ ^ \ , '-l^^ . (10) 



There arc two reasons why we prefer to study the preliminary data for the asymmetries weighted with a 
power of Ph± [7] instead of the final data for the asymmetries weighted without Ph± [8] , in spite of the caveat 
that preliminary data can be subject to changes due to refined data analyses. Firstly, in the parton model 
approximation only asymmetries weighted with an appropriate power of transverse momentum (e.g., Phj_ in 
the case of Sivers effect in SIDIS) allow a model independent disentanglement of transverse parton momenta 
in the target and in the fragmenting hadron [11]. An analysis of A^^'^'^~'^^^ would inevitably be biased by our 
prejudice concerning the distribution of transverse parton momenta in the target and in the fragmentation 
process, while the use of the asymmetry (9, 10) allows to avoid this problem elegantly. Secondly, below 
we will be interested in discussing the Sivers effect in DY pair production in the COMPASS experiment at 
considerably higher energies. It has been argued that asymmetries weighted without transverse momentum 
could be subject to strong dilution due to Sudakov effects, while this effect could be minimized by weighting 
the SSA by an appropriate power of transverse momentum [55]. 

Using the fact that the unpolarized distribution /f (x) and fragmentation Di(z) functions are well known 
and parameterized (see, e.g., [52, 53, 56]), one could try to extract directly information on the Sivers function 
using the so-called purity method which is being pursued by the HERMES Collaboration [54] . Instead, we 
choose a different strategy here and fit the HERMES data [7] for which we employ the ansatz 

2^/its!dis(^) ^ ~^/its!dis(^) = Ax^i^ - x)^ , 2;/^2^3|^jg(a;) = . (11) 

Several comments are in order. Firstly, we assume that the Sivers function for antiquarks can be neglected 
in comparison to the one for quarks. Secondly, we imposed the condition (7) derived from the large- A^^ limit 
[26]. Both assumptions are severe constraints, but ~ given the size of the experimental error bars [7] - they 
can be expected to hold with sufficient accuracy for our purposes. In fact, we shall see that the present 
data are compatible with the ansatz. Of course, one should expect that future precision data may demand 
to relax these constraints. At the present stage, however, these assumptions are very helpful to reduce the 
number of unknown quantities to only one, namely the Sivers u-quark distribution. We also neglect strange 
quark effects. 



^Note that throughout this paper we also neglect contributions from soft gluons [19, 20, 21]. 
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In order to illustrate to which extent the HERMES data allow 
to constrain the parameters in the ansatz (11) we performed two 
fitting procedures, one (I) with the parameter B = 1 fixed from the 
very beginning, and another one (II) where both parameters A and 
B were kept free. Using for fi{x) and D^^z) the parameterizations 
[52] (or [53] which yields a negligible difference) and [56] at = 
2.5 GeV^ we obtain the fits 
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Figure 3: Sivcrs function according to 
Eqs. (11, 12) as obtained from a fit to the 
HERMES data [7], see Figs. 4a-c. The un- 
polarized quark distributions xfl^(x) at = 
2.5 GeV'^, rescaled by the factor (— 1)/10, are 
shown for the sake of comparison. 



with a comparable x per data point of about 0.4. The fitting func- 
tions are of different shape, see Fig. 3, but they describe the HER- 
MES data [7] equally well, see Figs. 4a-c. The scale for the Sivers 
function in Eq. (12) corresponds to the average scale in the HER- 
MES experiment of (Q^) = 2.5 GeV^. We remark that the fits (12) 
are mainly constrained by the 7r+-data. Leaving the (and/or tt~ ) 
data out of the fit does not affect the numbers in (12) significantly. 

Thus we see that the experimental accuracy of the data does not allow one to constrain more sophisticated 
ansatze with more than two free parameters. Considering the discussion of the large- a; behaviour in the 
previous section, we have been guided to the ansatz (11). However, one should keep in mind that we use this 
ansatz only in the region x < 0.4 covered by HERMES, so the precise shape of f^{x) in the limit x — > 1 is 
of no relevance for us. 

Let us confront the results of our fit to the z-dependent data from [7] . Since the latter was not used to 
constrain the fit, the comparison in Figs. 4d-f can be viewed as a "cross check" of the fitting procedure. The 
expression for the asymmetry is given by Eq. (10) but with the average with respect to x instead of z. The 
shape of the SSA is dictated by the parameterization for Di(z) from Ref. [56]. The asymmetry is linearly 
rising with z for tt*^ (where Df is the same for all q = u, u, rf, d) and nearly so for 7r+ (where favoured flavour 
approximation works well), but it has a peculiar shape for tt^ (where 1/iVc-corrections to the Sivers function 
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Figure 4: (a,b,c) The azimuthal SSA A^IJ!'''* »'s;-hj./--iv ^ function of x. The preliminary data are from the HERMES 
experiment [7]. The curves are obtained from the large-N^ constrained fits I and H (denoted as in Fig. 3) of the Sivers function. 

(d,c,f) Ay'^'"*'' 't's)Ph±/Mfi ^ function of z, with the preliminary data from [7], and the theoretical curves from the fits I and 
H of the Sivers function. The z-dependent data were not used for the fit, and serve as a cross check of our results. 
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would have the most impact). Wc conchide that the ansatz (11) and the fits (12) are well compatible with 
the z-dependence of the data, see Figs. 4d-f. 

We observe that the obtained fit satisfies \f^ll^^^{x)\ < ^/f(a;), see Fig. 3. Multiplying Eq. (3) by 

Pt|/(2Mjv) and integrating it over transverse momenta gives the inequality \fiT]iii}isi^)\ — 2mI- fii^) 
which (defines a phenomenological mean parton transverse momentum {px) and) is less restrictive than the 
inequality observed in Fig. 3, if we assume (px) ~ 0.8 GeV (see Ref. [44]). In this sense, we note that our 
result is in agreement with the positivity bound in Eq. (3). 

Let us also remark that the HERMES data [7] are compatible with the large- A^c counting rule in Eq. (7) 
within their present statistical accuracy, which is proven by the fact that a fit with the ansatz (11) works. 
The sum rule (6) is satisfied by our parameterization (recall the suppression of the gluon Sivers function in 
the large- iVc limit). 

Experiments with the deuterium target, for which'' frr^^ ~ frr^^ + /ir*^" = frr + fir^ ^^c., are 
best suited to study deviations from the ansatz (11). Thus, the Sivers effect spin asymmetries from the 
deuterium target are suppressed with respect to proton asymmetries by a power of Nc in the large- A^c limit. 
The preliminary COMPASS data on A^^'^" '^^ '^ from a deuterium target are compatible with zero within 
error bars [9], and thus do not contradict the large- A^c motivated ansatz (11). 

Let us compare our result to parameterizations of the Sivers function in the literature obtained from 
studies of SSA in p^p ttX. The fits in Eq. (12) agree in sign, but are shifted towards smaller x, and one 
order of magnitude larger than the Sivers function obtained in Ref. [43]. On the basis of the latter it was 
estimated [40] that the Sivers effect can be neglected with respect to the Collins effect in the twist-3 SSA 
^uif observed at HERMES [4]. Our considerably more sizeable result does not support these conclusions 
and indicates a possible necessity to reconsider the interpretation [57, 58, 59] of the HERMES data on 
though recent studies indicate that the Sivers effect cannot play a dominant role in this SSA [60]. 

In the earlier extractions [12, 43] fir^^'^ix) was so small due to a small value of (pt) < 0.2 GeV assumed 
in these references. With the updated value {pr) ~ 0.8 GeV the resulting updated Sivers functions [44] 
which fits the FNAL data [2] is of comparable magnitude as our result. However, the different a;-shape of 
the /i^'s'iDis required to describe HERMES data as compared to the /{^^^^' required to describe the FNAL 
data remains. 

Finally we note that the prediction of the sign of the Sivers function for u- and d-flavour in Eq. (8) is 

confirmed by our parameterization, which suggests that the physical picture of the Sivers effect of Ref. [47] 
is apparently able to catch main features of the Sivers effect. 



4 Sivers effect in the Drell-Yan process 

The information on deduced in Sec. 3 from the HERMES data [7] is rough, however, as we shall see, 
sufficient for our goal to predict the sign and to gain insight into the magnitude of the SSA in DY. 

The process p^h fi'^fj~X (with h = p, 7r~ in the following) is characterized by the variables s = 
(pi +P2)'^, the dilepton invariant mass = {ki + k2f' with (and fci/2) indicating the momenta of the 
incoming proton and hadron h (and the outgoing lepton pair), and the rapidity 

1 , Pi{ki + fe) 

^=2^^^^(fc7TM- ^''^ 

Let us consider the azimuthal SSA which is weighted by |qT|, the dilepton momentum transverse with respect 
to the collision axis, and defined as a sum over the events i according to 



^sin(0-0s)^ ^ Ei ^H^i - <l>s,i)^-^{N\<i>u <Ps,i) - Ni-{cPi, cl>s,i + n)} 



where t (i) denote the transverse polarization of the proton. (See Fig. 2b for the definition of the kinematics.) 
To leading order the SSA is given by 

sin(,-,.)^ _ Eaelx,fl-^'^:/'{x,)x,f-^/\x,) 



•^Wc neglect nuclear binding effects and assume isospin symmetry which is legitimate given the present level of accuracy. 
Parton distributions without a target label (D = deuteron, p = proton, n = neutron) refer, as everywhere in this note, to the 
proton. 
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Figure 5: The azimuthal SSA A^^^'^'" 'I's)i±/Mn Drell-Yan lepton pair production, pT h fi+pi-X, as function of y: (a) 
for the kinematics of the PAX experiment where the hadron h = p, (b) for the kinematics of the COMPASS experiment where 
h = 7r~ . The different curves correspond to the fits I and II (see Eq. (12)), including the sign-reversal in (1). 



where the parton momenta Xi/2 in Eq. (15) are fixed in terms of s, Q'^ and y, 



Xl/2 = 




(16) 



The sums in Eq. (15) run over all quark and antiquark flavours, and we indicate explicitly to which hadron 
the distributions refer. 

In the PAX experiment antiprotons with a beam energy of 25 GeV could be available, i.e. s = 45 GeV^. 
In this kinematics one could explore the region aroimd = 2.5 GeV^, which is below the region of J/^ 
production, and well above the region of dileptons from $(1020)-decays. Taking into account the change of 
sign in the Sivers function in DY as compared to SIDIS, see Eq. (1), we obtain the result shown in Fig. 5a.* 
We observe that the two fits I and II, which describe the HERMES data of SIDIS equally well, give clearly 
distinguishable results in DY. Considering depolarization, detector acceptance and other effects, it might be 
difHcult to distinguish the effect of the different parameterizations in Eq. (12). However, the asymmetry is 
large enough to check unambiguously the QCD prediction of the different sign of the Sivers function in DY 
and SIDIS. 

In the COMPASS experiment using a 7r~ beam (s = 400 GeV^) one could also measure the asymmetry 
(15). In Fig. 5b we show the asymmetry for = 20GeV^ using for the pion the parameterization from 
Ref. [61]. Although ft^^ix) is far less constrained by data compared to /°^^(x) the result in Fig. 5b is 
rather insensitive to the choice of parameterization, and changes very little if we use the pion distributions of 
Ref. [62] (consistently in combination with the nucleon distributions from Ref. [53]). We observe a situation, 
which is qualitatively and quantitatively similar to the case of DY from pp-collisions. Note that we neglected 
evolution effects (from = 2.5 GeV^ in Eq. (12) to = 20GeV^ in Fig. 5) for the Sivers function. 
However, the influence of evolution is presumably much smaller than other uncertainties in our study. Note 
that by using the g^-weighted SSA we have avoided another serious problem in this context, namely Sudakov 
suppression [55], see the remarks in the previous section. 

In order to extract quantitative information from the future COMPASS and PAX experiments it is 
necessary to go beyond the LO formalism, to consider effects of soft gluons and K-factors, and to study the 
role of possible higher twist effects. The corrections due to these effects cannot be expected to be negligible. 
However, they are unlikely to be able to change the sign of the asymmetry. Thus, both the COMPASS as 
well as the PAX experiment could provide a thorough experimental test of the QCD prediction in Eq. (1). 

SSA in DY can also be studied at RHIC in p^p n'^ ii~ X. Since only one proton needs to be polarized 
the counting rates would be somehow more sizeable than in the case of double spin asymmetries related to 
the transversity distribution hi{x) which are, however, small [63]. Moreover, in this case, one is sensitive to 
the Sivers antiquark distribution which is not constrained by the HERMES data. We remark that the RHIC 
experiment is well suited to learn, e.g., about the Sivers function from SSA in p^p — > -kX [3] or the gluon 
Sivers function [64, 65]. 

^Thc DY asymmetry appears positive like the SIDIS asymmetry at HERMES, despite the change of sign of the Sivers 
function due to conventions: In DY we define the z-axis in the direction in which the polcirized particle moves. In SIDIS at 
HERMES the z-axis is defined in the opposite direction, see Figs. 2a and 2b. 
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5 Conclusions 



The recently reported HERMES data [7, 8] on SSA provide a theoretically unambiguous experimental ev- 
idence for the existence of T-odd distribution (and fragmentation) functions. We analyzed the HERMES 
data and demonstrated that they are consistent with predictions from the large- A^c limit of QCD [26] for the 
Sivers fiuictions, namely = modulo 1/Nc corrections. Imposing this large- A^,. result as an exact 

constraint we were able to obtain parameterizations of the Sivers quark distribution functions. The neglect 
oil/Nc corrections (as well as antiquark effects) in a first approximation is reasonable, keeping in mind the 
large error bars of the present data which do not allow to constrain more sophisticated ansatze. 

On the basis of the obtained parameterizations we estimated SSA in the Drell-Yan process for the PAX 
[p^p — > iX^^~X) and COMPASS {p^-K~ /z+/i~X) experiment. According to the theoretical understanding 
of T-odd parton distributions in QCD the Sivers function should obey a particular universality relation, 
namely appear with opposite sign in DY and SIDIS [17]. Our estimates show that both experiments could 
be able to test this prediction, which would be a crucial check of the present understanding of T-odd 
distribution functions and the QCD factorization approach to the description of SSA. 
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Note added: After this manuscript was completed the COMPASS collaboration has published final data 
on transverse target SSA from a deuteron target [66]. We stress that the large- prediction for the flavour 
dependence of the Sivers function [26] naturally explains the compatibility of the sizeable Sivers SSA from a 
proton target observed at HERMES [7, 8] and the small (consistent with zero within error bars) SSA from 
a deuteron target observed at COMPASS [9, 66]. The Sivers effect in the deuteron is sensitive to the flavour 
combination (/fy* -|- /ir*), and thus suppressed with respect to the Sivers effect in the proton by one order of 
Nc in the large- A^c expansion. In nature A^c = 3 is sufficiently large to explain the observations - considering 
the statistics of the first experiments [7, 8, 9, 66]. 
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